Andreev states near short-ranged pairing potential impurities 
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We study Andreev states near atomic scale modulations in the pairing potential in both s- and 
d-wave superconductors with short coherence lengths. For a moderate reduction of the local gap, the 
states exist only close to the gap edge. If one allows for local sign changes of the order parameter, 
however, resonances can occur at energies close to the Fermi level. The local density of states 
(LDOS) around such pairing potential defects strongly resembles the patterns observed by tunneling 
measurements around Zn impurities in Bi2Sr2CaCu2 08+a: (BSCCO). We discuss how this phase 
impurity model of the Zn LDOS pattern can be distinguished from other proposals experimentally. 

PACS numbers: 74.45.+C, 74.72.-h, 74.62.Dh 



Motivated by the experimental ability to determine the 
LDOS with high resolution in both energy and real space 
using scanning tunneling microscopy (STM), there has 
recently been a large interest in the perturbations caused 
by impurities in superconductors. This is because impu- 
rities disturb the underlying superconducting state and 
hence constitute a natural probe of the state in which 
they are embedded Ij. For magnetic adatoms on the sur- 
face of conventional s-wave superconductors such exper- 
iments were performed by Yazdani et aiQ ■ In the super- 
conducting state of BSCCO, STM measurements have 
provided detailed information about the electronic struc- 
ture near Ni and Zn impurities 0, 0, 1^. Near Zn it was 
found that each impurity generates a sharp conductance 
peak close to the Fermi level {lub ~ — l.SmeV). By fixing 
the bias voltage between the tip and the sample ai ujb, 
the spatial structure of this state can be mapped out: it 
is strongly localized near the Zn atom with maximum in- 
tensity on the impurity site and second largest intensity 
on the next-nearest neighbor sites. 

Although there exists a large number of theoretical 
treatments dealing with the resonant states generated 
by nonmagnetic impurities in d-wave superconductors, 
no consensus has been reached on their relevance to 
experiments 0. Within the most straightforward sce- 
nario where Zn is modelled as a delta-function potential 
scatterer, low-energy resonant states are generated in the 
unitary limit where the potential is large compared to all 
other energy scales of the problem 'Gj. For realistic bands 
with particle-hole asymmetry, a particular value of the 
impurity potential V must be chosen to tune the reso- 
nance energy lub to be close to the Fermi level. Thus, the 
potential scattering model can reproduce the energetics 
of the resonant state. However, as is well-known, it pro- 
duces a low-energy spatial LDOS pattern with a severely 
suppressed amplitude on the impurity site, in contrast 
to experimental results on Zn 5]. Physically, this is clear 
since a large on-site potential V penalizes any substan- 
tial amplitude of the impurity-state wavefunction on the 
impurity site. In fact, within this model the largest inten- 
sity is found on the nearest neighbor site to the impurity. 



These properties of the delta-function potential scatterer 
are shown in Fig^ The discrepancies between experi- 
ments and the potential scattering scenario have led to 
the suggestion of various filter functions motivated by the 
fact that the measurements are performed on the top- 
most BiO layer whereas the nonmagnetic impurity sub- 
stitutes a Cu atom in the Cu02 plane two layers below 
the top BiO laver|5 IgUlpf. The significance of the filter 
functions remains controversial. Other scenarios have fo- 
cused on the formation of local moments near the Zn 
impuritieslHElEl- 

For instance, if the main effect of 
nonmagnetic impurities in cuprate superconductors is to 
locally break spin singlet bonds, the sharp conductance 
peak has been interpreted as a Kondo resonance |l2j. 

While most models assume that an impurity produces 
a screened electrostatic potential and possible magnetic 
exchange effects, little attention has been paid to the 
possibility that a defect can cause local distortions of the 
interactions which lead to BCS pairing. We proposed 
this idea recently in a phenomenological approach to un- 
derstanding the correlations between dopant atoms and 
gap inhomogeneity in Bscco[i3. In that paper, it was 
shown that by assuming that out of plane dopant atoms 
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FIG. 1: (Color online) Left: LDOS in a d-wave supercon- 
ductor at the impurity site for varying potential scattering 
strengths V (clean case: dashed curve). As V increases the 
resonant state sharpens and moves to lower energy. Right: 
Real-space LDOS pattern at the resonant energy ujb with a 
delta- function impurity potential {V = 5t) at the centerQ. 
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locally enhance the pair interaction, many of the peculiar 
experimental correlations could be understood (2^. Here, 
we continue this phenomenological approach, and calcu- 
late the Andreev states associated with atomic scale per- 
turbations in the pairing potential relevant for supercon- 
ductors with very short coherence lengths. This is con- 
trary to the conventional study of Andreev states within 
the quasiclassical approach where physical quantities are 
assumed to vary slowly on the atomic scale. From the 
results for the d-wave case, we show that the low-energy 
conductance peaks observed near in-plane nonmagnetic 
impurities in BSCCO can be explained in terms of local- 
ized Andreev states. We further investigate the robust- 
ness of this picture and suggest some measurements that 
could be used to test this scenario. 

The model is given by the usual BCS Hamiltonian 



Wo = ^ SkcL5k. + (^kc[^cLki + H.c 



(1) 



k.o- 



where = — 2t(cosfca; -I- cosfcj^) — 4t' cos /c^: cos fcj^ — /i 
denotes the quasiparticle dispersion with nearest (next- 
nearest) neighbor hopping t(t'). For s-wave pair- 
ing, Ak — Aq, whereas for d-wave pairing, Ak — 
^{coskx — cosfcj,). In terms of the Nambu spinor 

ipl^ = (c]jj.|, c_ki), the corresponding Green's function 
uJnT)) in Matsubara representation is given by 

G [K ZW„) = -2 , (2 

where — ^k + ^ki ^^^d denote the Pauli matrices. In 
real-space, the perturbation due to delta-function poten- 
tial impurities in the diagonal channel and modulated 
pairing in the off-diagonal ri channel is given by 



n'{r, r') = [VS{r)Sir')r3 ~ 5A{r, r')ri] 



(3) 



To obtain the resulting LDOS as a function of energy 
and lattice sites, one needs to determine the full Green's 
function CJ(r, icun) given by the Dyson equation 

Gir, r') = ty°(r - r') + ^ Cy(r, r")i?'(r", r"')t/°(r"' - r'). 

(4) 



Thus, by calculating the matrix elements of 



g°(r,.....)^i: (-"-^+^rtf'--^) exp(.k.r), (5) 



k 



the remaining problem is that of a matrix inversion. The 
solution is presented in terms of the so-called T-matrix 

giv,v') = g°(r-r') + ^g°(r-r")T(r",r"')^°(r"'-r'). 

(6) 

The poles of the T-matrix, or equivalently, the determi- 
nant of (1 — H'G'^), determine the bound state energies. 



For s-wave superconductors, it is well-known that non- 
magnetic impurities cannot generate states inside the 
gap[l3- However, as shown by Yu and Shiba, point-like 
magnetic impurities {H'{r,r') = i/)! [Kni5(r)(5(r')To] i/^r') 
of strength Vm produce bound states at ujb = ±Ao(l — 
(Ki7i-7V(0))^)/(1 + {V„,T:N{0)f), where N{0) is the den- 
sity of states at the Fermi level in the normal state 17, 1^. 
Assuming, as a toy-model, a point-like ri scatterer in 
an s-wave superconductor, a straightforward calcula- 
tion shows that an attractive delta-function ri scatterer 
generates Andreev bound states with energies given by 
LUB = ±Ao|(l - (7r7V(0)<5A)2)/(l + {TrN{0)6A)^)\. Here, 
6 A is treated as a free parameter and not necessarily just 
the self-consistent response of the order parameter to e.g. 
an electrostatic potential. A cursory examination of this 
expression for the bound state energies ub shows that 
as the local gap decreases from the bulk value Aq, the 
Andreev states are pushed further into the gap. While 
initially located near the gap edge, for increased scatter- 
ing strength the bound states can in principle be tuned 
all the way to zero energy, but this requires SA of order 
the Fermi energy, outside the framework of conventional 
weak-coupling BCS theory. 

What about d-wave superconductors? For a single 
nonmagnetic T3 delta-function impurity, the T-matrix 
is diagonal with Tn = S{r)[V-^ - ^LH^]-^{r'), 

T22 = <5(r)[-y-i - Ek 7j^]"'^(i-')- It is the poles of 

k 

these matrix elements that largely determine the proper- 
ties presented in Fig^ We expect short-ranged ti scat- 
terers to be more relevant in the case of d-wave cuprate 
superconductors due to the short coherence length of 
these materials. In the d-wave case, a point-like ti im- 
purity is defined as a modulation of the pairing potential 
on the four links attached to the impurity site. Following 
Shnirnian et al. 15], in the case of a particle- hole symmet- 
ric band, one can determine the subgap poles from the 
determinant D{uj) given by 

D{uj) = 1 - aL{uj) (2 - aL{uj) + aujP{uj)) , (7) 

where a parameterizes the local gap suppression, 6 A — 
aAo. Here, P{lu) and L{uj) are given by 



(PM,LH) = - 



k 



f2 

Sk 



(8) 



We can rewrite Eq. (0 in terms of: P{u!) = loK{1/{1 — 
Cu^ )/{VT^ ), L{uj) = Cuy{VT^)K{l/{l - Cu^)) + 
(Vl ~Cj^)E{l/{l-Cj^)), where E{x) {K{x)) is the com- 
plete elliptic integral (of the first kind), and w = lo/Aq 
and a — 47rA^(0)Aoa. In Fig|2^ we plot the real part 
of the determinant, Rel?(a;), for different values of the 
strength of the off-diagonal impurity potential. The 
imaginary part, 1tiiD{uj) (not shown), is a monotoni- 
cally decreasing function as ^ with ImZ?(0) = 0. 
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FIG. 2: (Color online) (a) Real part of the determinant, 
ReD{Lu), with a = 0.4 (top), a = 0.6 (middle), and a = 0.8 
(bottom), (b) LDOS at the impurity site for the same values 
of a as in (a), the dashed curve displays the clean result. 

Thus, the point-like ri impurity in a d-wave supercon- 
ductor also supports Andreev resonant states. The re- 
sulting LDOS at the impurity site is shown in Fig|23: for 
increased strength of the gap suppression the Andreev 
resonance sharpens up and moves to lower energy. As 
opposed to the T3 point-like impurity, the resonance ex- 
ists symmetrically around zero energy with strongly in- 
creased amplitude as lob — > 0. In general, this symmetry 
is broken by including e.g. weak T3 potentials. One ex- 
pects the same to be true for more realistic non-particle- 
hole symmetric band structures. In this case, however, it 
is advantageous to resort to numerics. Since we will be 
mostly interested in the case of optimally doped BSCCO, 
in the following we fix t' = —0.3t, Aq = OAt, and ^ — —t 
corresponding to 17% hole doping. 

In FigOlwe show the LDOS for this band structure at 
the three sites nearest to the impurity. Clearly, in this 
case the particle-hole symmetry is broken, and the nega- 
tive bias resonance dominates the LDOS at the impurity 
site. Figl^^d) shows explicitly that the ti scattering chan- 
nel produces the largest maximum on the impurity site 
and has a second maximum on the next-nearest neighbor 
site. If we apply the filter function of Ref. 0, the LDOS 
pattern becomes similar to the one shown in Fig^Ib) in 
disagreement with experiments. When including the T3 
scattering channel the asymmetry in FigsOKibjC becomes 
more pronounced. Depending on the specific band struc- 
ture, we find that in general the main features of Fig|3 
remain valid for V t. Therefore, if the conventional 
argument is eventually proven correct that a closed un- 
screened d-shell on Zn necessarily implies a local poten- 
tial of several eV, then the present model cannot be ap- 
plied to explain the low-energy conductance peak seen 
near these impurities in BSCCO 0. 

The important point of the results presented here is 
that point-like ti scatterers generate well-defined An- 
dreev states in both s- and d-wave superconductors. As a 
specific example (FigEJ, we found that in order to model 
the sharp zero bias conductance peak near Zn impurities 
in BSCCO, the local gap has reversed sign compared to 
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FIG. 3: (Color online) LDOS for various strengths of the ri 
scatterer at the impurity site (a), on the nearest neighbor site 
(b), and the next-nearest neighbor site (c). In each figure, 
the dashed curve shows the bulk LDOS. (d) The real-space 
LDOS distribution at the resonance energy ujb when a = L4. 



its value in the bulk[l9j. Thus, it is intriguing to consider 
the idea that nonmagnetic planar impurities in cuprate 
superconductors act mainly as phase impurities changing 
locally the sign of A on the x and y links compared to 
the bulk. The particular value SAb necessary to generate 
well-defined low-energy Andreev states near the impurity 
site depends on band structure, possible nonzero T3 scat- 
tering, and the spatial extent of the pairing modulation. 
Regarding the last, we find that for longer ranged phase 
impurities, the necessary strength SAb is substantially 
reduced. For example, for a ti impurity ranging over 
two lattice constants, a = 0.5 generates LDOS similar to 
the solid black curves of Fig|31 Thus, 6Ab of order Aq 
is sufficient to generate well-defined low-energy states. 

Though the approach here is phenomenological and 
the microscopic origin of these phase impurities is un- 
known, they clearly must be generated by the local per- 
turbations resulting from a Zn ion replacing a Cu. In 
recent work, it was found that the main features of the 
LDOS spectra in the inhomogeneous nano-scale regions 
of BSCCO materials can be accounted for by local dopant 
induced modulations of the pairing interaction[lj,l20j. In 
the present context, we find that self-consistent solutions 
of the Bogoliubov-de Gennes equations generate phase 
impurities if the pair interaction near the Zn impurity 
position is repulsive. Local pairing modulations can be 
generated, for instance, by locally modified magnetic ex- 
changes and/or electron-phonon couplings. A similar ap- 
proach has been used previously to explain the enhanced 
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FIG. 4: (Color online) LDOS at (0, 0) (a) and (2, 0) (b) when 
two impurities are positioned at ( — 1, 1) and (1, —1) for both 
the Ta (solid) and ri (dashed) channel. In both cases, the 
T3 LDOS has been rescaled for clarity, (c-d) LDOS at the 
impurity site for a single T3 (c) and n (d) scatterer of fixed 
strength (V = 5t for rs, and a — lA for n) but varying Ao. 



Tc in twinning-plane superconductors (2 Ij . 

We end this section by suggesting alternative experi- 
mental tests of whether the ti or T3 channel is dominating 
the LDOS near Zn in the BSCCO materials. In both sce- 
narios the appearance of the resonance peak is tied to the 
superconducting transition temperature Tc. Recently, 
there has been an increased interest in the quantum in- 
terference caused by nearby impurity states |2a. I23L l24l|. 
Figs0jL,b show the resulting LDOS from two interfering 
T3 or Ti impurities. Here, the two impurities are fixed in 
close proximity at (—1,1) and (1,-1) and the LDOS is 
shown at (0,0) (a) and (2,0) (b). Evidently, there is a 
qualitative difference in the resulting interference pattern 
between the two scenarios: whereas the T3 resonances are 
absent, the well-defined Andreev resonances are split and 
enhanced by the interference. The latter point remains 
valid for other sites (except at the impurity site) in the 
vicinity of the two Andreev resonances. Another possi- 
ble experimental test utilizes the ubiquitous gap inhomo- 
geneities observed near the surface of BSCCO |25lERl27j . 
Specifically, as shown in Figs01:,d the LDOS near Zn 
impurities in large and small gap regions is different in 
the two scenarios: whereas scatterers change both lub 
and the resonance amplitude, the ti case exhibits mainly 
amplitude modulations with the largest resonance ampli- 
tude for intermediate gap Aq. 

In summary, we have studied the states near short- 
ranged off-diagonal impurities in both s- and d-wave su- 
perconductors. In both cases we find that the impurities 
generate low-energy Andreev states. This offers a new 
possibility for the origin of the low-energy conductance 
peak observed near Zn impurities in BSCCO . 
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